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RADIATOR DESIGN AO INSTALLATION 



By M. J. Brevoort and M. Leifer 



SUMMARY 



The fundamental principles of fluid flow, pressure 
losses, and heat transfer have "been presented and analyzed 
for the case of a smooth tube with fully developed turbulent 
flow. These equations apply to tubes with large length- 
diameter ratios where the flow is at a high Reynolds Number. 
The error introduced "by using these equations increases as 
the magnitude of the tube length and the air-flow Reynolds 
Number approaches the values encountered in modern radiator 
designs. Accordingly, heat-transfer tests on radiator sec- 
tions were made and the results are presented in nondimen- 
sional form to facilitate their use and for comparison with 
other heat-transfer data. In addition, pressure losses were 
measured along smooth tubes of circular, square, and rectan- 
gular cross section and the results were also correlated and 
are presented in nondimensional form. 

The problem of a radiator design for a particular in- 
stallation is solved,, the experimental heat-transfer and 
pressure-loss data being used, on a basis of power chargeable 
to the radiator for form drag, for propelling the weight, 
and for forcing the air through the radiator. 

The case of an installation within a wing or an engine 
nacelle is considered. An illustration of radiator design 
is carried through for an arbitrary set of conditions. Suf- 
ficient detail is given to enable the reader to reproduce 
the analysis for any given case. The effect upon the cool- 
ing efficiency is considered for a change in the following 
design conditions: 



(1 
(2 

(3 
(4 

(5 
(6 
(7 



The pressure drop across the radiator. 

The pump efficiency of the duct. 

The free-area ratio of the radiator. 

The wing characteristics for the design condi- 
tions. 

The air constants. 

The change in density with altitude. 

The width of the water passageway. 
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Comparisons of the cooling efficiency possible with 
tubes of various shapes lead to the conclusion that the 
smooth hexagonally shaped tube is superior to the circular 
tube and that finned tubes and tubes shaped to produce 
turbulence cause a decrease in cooling efficiency. 

The use of ethylene glycol instead of water for the 
liquid makes it possible to increase the cooling efficien- 
cy several times (depending upon the inlet-air tempera- 
ture). The analysis of the length of passage on the 
liquid side of the radiator shows it to be a secondary 
consideration. The design problem for an installation in 
a separate nacelle is solved. Several values are chosen 
for the effective nacelle-drag coefficient to permit evalu 
ation of the nacelle-wing interference; on the cooling effi- 
ciency. A cooling-power comparison between liquid-cooled 
and air-cooled engines is made and it is concluded that 
the power to cool is not a determining factor in the selec* 
tion of the type of engine. The recovery of mechanical 
energy from the heat energy dissipated by the radiator in 
the duct is shown to be possible. An approximate deter- 
mination of the magnitude of the effect is made. The de- 
sign of the radiator is not affected by considerations of 
the energy recovery. 

INTRODUCTION 



With the wealth of information available in heat- 
transfer processes, all the necessary information for the 
design of aircraft radiators would appear to be at hand. 
Such is practically the case. In fact, there is so much 
information on special cases that a designer setting out 
to choose a radiator is confused by a large variety of re- 
ports, all. of which appear to have some bearing on the 
problem. 

The most misleading reports of this type are those 
dealing with, the radiator mounted in the air stream. The 
data in these reports were invaluable when radiators were 
mounted outside the airplane without benefit of cowling, 
but the data apply only to a special case and cannot be 
used for a modern radiator installation in which the radi- 
ator is located in a wing, a cowling, «r an engine nacelle 

A radiator mounted in a free air stream acts more or 
less like a flat plate, depending upon the resistance to 
flow through the radiator tubes. If the radiator has a 
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high resistance, most of the air flow approaching the radi- 
ator diverges and passes around the outside. In so doing, 
the air "breaks away at "both the front and the rear faces 
of the radiator, resulting in a large drag. If the radia- 
tor is made with a low resistance to flow through the tubes, 
most of the air approaching the radiator goes through the 
radiator and the drag is greatly reduced. 

When the radiator is mounted in the free air stream, 
the power consumed "by the radiator increases with the cube 
of the velocity. With the radiator in a duct, however, it 
is necessary to consider only the weight of the radiator 
and the internal power consumed in it. The radiator locat- 
ed in a duct can he made to dissipate a given amount of 
heat "by almost as low an expenditure of power for cooling 
as desired. But the lowor the power, the larger the radia- 
tor. The most economical radiator size dopends upon the 
wing loading of the particular airplane. 

There appears to he a prevalent notion that, "by the 
invention of some peculiar tube arrangement and shape, 
large gains in cooling performance can he attained. The 
"best situation that can he obtained is to have the fric- 
tion loss that occurs in the radiator he true surface fric- 
tion on the direct cooling surface. Pressure losses asso- 
ciated with sharp contractions, expansions, and sharp hends 
as well as surface friction over an indirect cooling sur- 
face are all undesirahle. If the optimum condition is re- 
alized, it appears to ho of little importance from consid- 
erations of hoat dissipation whether the tuhes are circu- 
lar, square, rectangular, hexagonal, or any other conven- 
ient shape. One shape may have an advantage over another 
as regards the ratio of weight required to the space re- 
maining for the liquid flow. If the radiator meets the 
requirements of direct cooling surface and pressure loss 
connected only with the surface friction, however, the 
chief consideration may he the ease of manufacture, the 
ease of repair, and the durability. 

This report will review the fundamental principles of 
heat transfer as they apply to the problem of radiators. 
Original data on radiators will be presented and a number 
of special cases will be analyzed on the basis of the fun- 
damental principles, the data, and similar data from, other 
source s . 
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SYMBOLS 



The following symbols are used in the report and are 
listed alphabetically for ready reference. 



A 



P 
C V 



'D 



"a 

H 

H 



frontal area of the housing. 

total cooling surface upon which h^. is based, 
constant s . 

specific heat at constant pressure, 
specific heat at constant volume, 
lift coefficient of the wing, 
drag coefficient of the wing. 

drag coefficient of the radiator installation 
with no cooling-air flow. 

hydraulic diameter of the tube. 

drag. 

radiator depth. 

/Ap D x 

friction factor ( ~— ) 

free-area ratio, of open frontal area to the 
total frontal area. 

the acceleration of gravity. 

surface heat-transfer coefficient. 

over-all heat-transfer coefficient from fluid 
to fluid. 

heat-transfer coefficient from air to tube wall. 

heat-transfer coefficient from liquid to tube 
wall. 

quantity of heat dissipated per unit time. 
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k, 
M, 



Nu = 



fe2 

k 



P 

Ap, 

Ap f . 

Pr=^ 

k 

q. =4 p V 

Q 



a. 
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T 



la 

T 



required quantity of. heat dissipation, 
thermal conductivity, 
constants . 
tube length. 

mass flow of fluid per unit time. 

Husselt number. 

power. 

total power chargeable to the radiator. 

power required to force the air through the 
radiator . 

power required to support and propel the weight 
of the radiator. 

static pressure. 

total pressure difference, usually across the 
radiator . 

drop in pressure due to skin friction. 
Prandtl number . 

dynamic pressure (q carries same subscripts as 
V, velocity). 

quantity of air by volume passing through the radi- 
ator per unit' time. 

Reynolds Slumber. 

universal gas constant. 

absolute temperature, °P. 

air temperature. 

inlet-air temperature, 
average liquid temperature. 



Aa? a = T a - Tj_ a , change in air temperature. 

AT W , change in liquid temperature. 

^ T w-a' mea:a temperature difference "between coolant 
and air. 

V^, air velocity in the duct. 

7^, air velocity in the tube. 

Vq , air-stream velocity. 

v, specific volume. 

W, the work obtainable from a cycle. 

w, radiator width. 

W r , weight of the radiator. 

x, length of tube to the section under consider- 
ation. 

e, factor by which to multiply radiator weight 
to account for the additional required 
airplane structure. 

V, ratio of the specific heats of air. 

p, density, 

p^, density at altitude, h. 

v, kinematic viscosity. 

T) t , heat-transfer efficiency. 

"Hp, 'pump efficiency of the duct with the radi- 
ator installed. 

coefficient of viscosity. 

Subscripts : 

i, apply to the inside of the tube. 
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o, apply to the outside of the tube. ! . 

0, apply to the free air stream. 

t, apply to the liquid side of the radiator. 



S INGLE TUBES 
Analysis 



A radiator "being composed of a number of small tubes, 
it is essential to have a clear understanding of the op- ■ 
eration of the individual tube in a study of radiator 
characteristics. This statement does not mean that single 
tubes' ©an be studied and the results applied without ' qual- 
ification to a radiator comprising an assembly of tubes. 

Huch of the published data on single tubes, has been 
taken for a region of the tube well downstream from the ■• 
tube entrance and, as a result, gives results that show 
the frictional losses and the heat transfer for tubes in 
which the boundary layer is in equilibrium. That is, the 
boundary layer is either in the laminar condition with ■ . . 
parabolic velocity distribution over the tube or is in the 
completely developed turbulent condition.. Now in a radi- 
ator tube, if the Reynolds Number is below some critical 

jal/ae., the boundary layer will be laminar and the velocity 
distribution will go through a continuous evolution from 
almost uniform at the entrance to parabolic across the 
section considered. If the tube is less than 150 diameters 
in length, the tube may be completely in the stabilizing 
region. It will therefore have a thinner boundary layer 
and, as a result, a greater friction and consequently a 
greater heat-transfer coefficient. If the Reynol ds, Nu mb er 
is above some cri tical value . the flow in the tube may be 
turbulent. If the entrance is favorable for initiating 
turbulence in the flow, the tube, will have turbulent flow 
throughout its length and will have the friction and the 
heat-transfer coefficients associated with turbulent flow 
and the particular Reynolds Number. If the entrance is 
designed to give a smooth streamline flow, the flow will 
evolve through a transition and become turbulent down- 
stream. The higher the Reynolds Number, the nearer this 
transition point will be to the -tube entrance. Thus it is 
seen that, even though the flow, may be above the critical 
Reynolds Number, the tube may be too short to c;ive an op- 
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portunity for the stabilizing flow to break down. If the 
entrance is such that turbulence is initiated, for a flow- 
above the critical Reynolds Number, the flow will remain 
turbulent throughout. On the other hand, if the flow is 
below the critical Reynolds Number, the turbulence set' up 
by an 'entrance will be damped out and laminar flow will 
re suit . 

It is now obvious that tests on single tubes hare di- 
rect application to the radiator problem only when the 
tubes are tested under conditions reproducing those ex- 
isting in a radiator. If the entrance to a tube assembly 
in a radiator is smoothly streamlined, then the heat- 
transfer coefficient decreases rapidly downstream from the 
entrance as the boundary layer builds up. It is there- 
fore necessary to study numerous combinations of lengths 
and diameters to develop the complete picture of the phe- 
nomena involved. If the entrance promotes turbulence, 
however, the friction and the heat-transfer coefficients 
will closely approximate the acceptod results for turbu- 
lent flow' in tubes. 

The pressure loss", or friction, in tubes has been 
studied by Stanton and others for both the laminar and the 
turbulent regimes in the regions where the boundary layer 
is developed. The relations that satisfy each condition 
are ^iven by G-lauert (reference l) as follows: 

-32 R~ 1 p —£~ for laminar flow (l) 



1 V+- 

-0.157 R~ % p— for turbulent flow (2) 



Friction-Loss Tests and Correlation 

The pressure losses were measured'on circular tubes 
having diameters of l/2, l/4, l/8, and l/l6 inch; a square 
tube of l/2 inch hydraulic diameter; and rectangular tubes 
1/2 by 1/8 inch and 1/2 by l/l6 inch. The length of tube 
was made appropriate to give a length-diameter ratio of 
150 to 200 for each tube. 

The pressures were measured by static tubes flush 
with the inside of the tube, as seen in figure l(a). The 



dp • 
dL 

dP 
dL 
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tube surfaces were smooth, and a hell— shape entrance as- 
sured streamline flow at the outset (fig. 1(h)). The 
Telocity was measured by a survey using a pitot tube mount- 
ed on a micrometer. Typical surreys for the l/2-inch- 
diaxoeter circular tube are shown in figure 2. It was 
found that the entrance was sufficiently close to a stream- 
line shape that the volocity calculated from the first 
static-tube reading and the atmospheric pressure agreed 
within the experimental error with the average velocity 
measured in the survey. 

Figure 3 shows a sample curve of the pressure-drop 
measurements on one of the tubes. The results for all the 
tubes tested, plotted in nondimensional form, are shown 
in figure 4; the abscissa is the Reynolds Number, the or- 
dinate is the friction factor f x , and each curve is for 

a constant length-diameter ratio. The data for all the 
tubes fall on the curves with the exception of the square 
and the rectangular tubes, which agree in the turbulent 
region but deviate radically from the data for circular 
tubes in the laminar region. 

From equations (l) and (2), f x = 16/R for laminar 
flow and .f a = 0.0785/R 0,85 for turbulent flow. The 
curves for fully developed laminar and turbulent flow are 
also shown on figure 4. The test data should fall between 
the two curves for fully developed flow because the sta- 
bilization and the transition regions are included. 

The friction-factor data for the radiators tested are 
shown in figure 5. Data from the single-tube tests at the 
same L/D ratio as that for the smooth-tube radiators are 
also, plotted. The curves from radiator data should fall 
higher than the curves from the single-tube data because 
the radiator Ap includes the losses at the exit. Since 
the transition occurs at the same Reynolds Number for both 
tubes and radiators, the exit losses can be determined 
from the differences between the friction factors for the 
single tubes and the radiators. ' Table I illustrates values 
of the exit l.oss over a range of Reynolds Numbers. The av- 
erage exit-loss value for the l/4-inch-diameter tube is 
seen to be- 0.2 , in agreement with Hartshorn (reference 

2). The larger values found for the l/8-inch-diameter 
tube may be attributed to a lower free-area ratio, which 
results in a greater expansion at the exit, and possibly 
to the relat-ively greater soldering irregularities at en- 
trance and exit. 
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TABLE I 

EXIT LOSSES FOR THREE RADIATORS 



Radiator 


Reynolds Number 


Exit loss/q.^ 




6 ,000 


0.19 


a A 


15 ,000 


.23 


30,000 


.305 




6,000 


0.18 


B 


15,000 


.15 




25 ,000 


.12 




6,000 


0.33 


C 


13,000 


.435 



See fig. 5 for explanation of the radiator designations. 



A radiator was made up having tube exits faired grad- 
ually from a circular to a hexagonal form over a distance 
of about 1 inch. Although this fairing appreciably reduced 
the exit loss for the single tube, no measurable improve- 
ment was found for the radiator. 

Figure 4 may be used to obtain the variation of air 
velocity in the tube with tube length for any pressure 
drop and tube diameter. The method is as follows: For 
any choice of I and D, figure 4 gives the relationship 
between f x = ApD/q^4L and Reynolds Number = p Y^.D/|j,. 

The conditions of temperature and pressure determine p 
and p.. The velocity V t is calculated from any Reynolds 
Number and follows. The corresponding value of f x 

gives Ap. Proceeding in this manner gives several sets 
of values of Y t and Ap. Then Y t is plotted against 

(Ap + 0.2 q t ) . The term 0.2 is included to account 

for pressure losses at the radiator exit. (See table I.) 
Such plots enable the dotermination of for any over- 

all pressure drop (= Ap + 0.2 q.^ ) • Repeating the process 
for other lengths or diameters will give the desired rela- 
tionship between air velocity in the tube and tube length. 
Figure 6 illustrates this relationship, using the values 
for D and Ap (over-all) shown on the figure. 
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The Reynolds Number describes the nature of the "bound- 
ary layer. Inasmuch as both the friction loss of the fluid 
and the heat transfer from wall to fluid depend upon the 
nature of the boundary layer, both characteristics are 
functions of the Reynolds Number. Reynolds recognized 
this common dependence on Reynolds lumber and stated the 
dependence in a principle known as Reynolds analogy. As a 
result 



Both the friction factor and the heat-transfer coef- 
ficient are larger with turbulent flow than with laminar 
flow above the critical Reynolds Number and each is larger 
with a thin than with a thick boundary layer. A more elab- 
orate relation has been evolved as a result of the work of 
Stanton, Prandtl, Taylor, and others that correlates the 
data on heat transfer from pipes for various velocities and 
fluids; namely, 

Hu = f = c, (Pr) 1 (R) m (§7 (4) 

Equation (4) is the equation used by many workers in • 
this field to correlate experimental heat-transfer data. 
The constant c and the exponents I, m, and n vary 
with a change in the type of -flow. For very ldng pipes, 
(L/D) n can be assimilated into the constant and a new 
constant c 2 used. The exponent n in any case is small. 

The term L/D has a strong influence only in the entrance 
region of the tube where the boundary layer is going 
through change. The Prandtl number is nearly constant 
whon air is used as the fluid flowing through the tube. 
It follows then that equation (4) for air flow and a con- 
stant L/D reduces to 

R m k ?t m P M , x 

h = c„ k or c 2 — ~" 5) 

a D 2 D i-m p_m 



Equations for Heat-Transfer and Cooling-Power Loss 

The local heat-transfer coefficient and the friction 
factor have been related in the foregoing discussion. 
There remains the problem of determining how much heat can 
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"be dissipated from a given tube as a function of wall tem- 
perature, inlet-air temperature, Reynolds Number, and 
length. Although the average local heat-transfer coeffi- 
cient can be determined by means of the tube friction, the 
tube length has an important bearing on the actual heat 
transfer by controlling the rise in the air temperature. 
The heat-transfer efficiency T| t has been defined to de- 
scribe the efficiency of heating up of the air as 

A -a 

n t = 

^w i ia 



It is desirable to have a high heat-transfer efficiency in 
order to dissipate a maximum amount of heat for a given 
mass flow of air. 

The heat dissipated and the power for cooling per 
unit open frontal area of tube can be determined as a 
function of Ap , T w - 2 i a > L » ^» an ^ * iie physical con- 
stants of the fluid. The ratio of these two quantities 
gives a relation from which it will be possible to deter- 
mine the optimum dimensions. 1 

The analysis will bo made for the case of funs'- de- 
veloped turbulent flow in a tube of constant wall tompor- 
aturc. The tube considered is circular, but the equations 
may be applied equally well to hexagonal or square tubes 
by using the hydraulic diameter. 

The heat AH dissipated for a length dx of the 
tube (fig. 7) is 



AH = h (T w - T a ) ttD dx (6) 



Also 



AH = p 7 t c p njfc. dTa (7) 
From equations (6) and (7) it may be shown that 



*a - 



■"•w na 



_ _4h x 

p T t c p D 

1& = T). = 1 ~ e (8) 



1 

This analysis for optimum tube dimensions takes no account 
of the radiator weight or the profile drag, which will 
be considered later. 
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where the heat-transfer -efficiency is now expressed as a 
function of the tube dimensions, the average heat-transfer 
coefficient, and the mass flow through the tube . 

For fully developed turbulent flow, Prandtl , quoted 
in reference 3 (p. 173), gives 



.? 53 °i c p Jots (P V t > 



(9) 



where a x = 0.0269. From equations (7), (8), and (9), a 
relation similar to that derived "by Hartshorn (reference 
2) is found: 



H = p T t c p 

where H 

P ,. 



ttD' 



< T w - T ia) 



1 - e 



o . a 



L 
D 



(10) 



is given in B.t.u./hr. 
in lb./cu.ft. 
in ft./hr. 

in B . t .u. /hr . / °F. / sq.ft. 
in lb./hr. ft. 



Temperatures are in degrees Fahrenheit and lengths are in 
feet. Hartshorn also developed a very useful family of 
curves of H per unit frontal area against pT at vari- 
ous values of I</D. 



From equation (2), the velocity is given as 



A 0 .571 -0.714 

y. = ±2 5 (11) 

t 0.571 0.57 1 0.14 3 0.571 

L c 2 T) p 

where c 2 = 0.157. 

Eliminating from equation (10), using equation 

(11), and taking the units into account, the heat dissipa- 
tion per unit open frontal area is given as: 
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H 



ttD / 4 



« ^, ~ 5 0.429' ■ 0.571 „0.714 , . 

0.841 x 10 c p p D . (T w - T ia ) 



/ 



x\l- 



c 0.571 ^0.143 ^0.571 



\ 



c c 0.1143 m0.3 pO.0285 £1,1143 



2.42 p 



0 . 086 



A 0.1142 .1.343 
&p D 



/ 



(12) 



where 



V 



is in ft . 3 /hr. 
Ap, in Id. /sq.ft. 



Figure 8 shows the variation of the heat dissipated 
with the length of the tube for a constant Ap for sev- 
eral tube diameters. Figure 9 shows the variation of the 
heat dissipated with the length of the tune for a tube 
diameter of l/48 foot for several values of Ap. It may 
he noted from figure 8 that the peak heat dissipation per 
square foot of open area occurs at about the same value of 
L/D for the various diameters. 



Hartshorn (reference 2) showed that the peak efficien- 
cy as a function of L/D varied as D 0 ' 2 . 



The basis for the choice of a radiator is its cooling 
efficiency, defined as the ratio of the heat dissipated to 
the- power chargeable to the radiator. 

The total power chargeable to the radiator is 

Ot) 0 At> 

P t = C D n A 1o \ + e cf. Y o W r + ^ 

On the right side of this equation", the first term is the 
power to overcome the form drag; the second term is the 
power to support the weight; and the third term is the 
power to force the air through the radiator. 



The pump efficiency is defined- in reference 4 for a 
separate nacelle installation as 
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The useful cooling work done per second is QAp and the 
work expended for cooling is (D - D 0 ) 7 0 . The quantity 
D 0 is the drag of a closed cowling with major dimensions 
similar to those of the actual co\vling, and (D - D 0 ) is 
the drag increment due to the cooling-air flow. 

Similarly, the pump efficiency for a duct is defined 
as the ratio of the useful cooling work done to the total 
work expended in forcing the air through the radiator. 
Thus for a wing duct, 

TL = 

v AD V 0 

where AD is the difference in drag between the original 
wing and the wing with the duct-radiator combination. In 
any installation, the entrance and the exit losses are in- 
cluded in the pump efficiency. 

The form-drag power is not entirely independent of 
the cooling-air flow, and the choice of the pump efficien- 
cy as an interference factor enables the . breakdown of the 
total power required into the subdivision stated. In an 
installation where the cooling-air flow may decrease the 
form drag, the pump efficiency may be more than 100 per- 
cent. It is possible that, for such cases, the basic shape 
without _ cooling-air flow may be considerably improved and 
such considerations may be worth investigating. 

It is likely that the pump efficiency will vary with 
the flight speed and the quantity of cooling-air flow. 

For the radiator mounted in a separate nacelle, the 
results presented in reference 4 give Cj^ and T|p ; the 

airplane design determines and C^. 

For the radiator mounted inside an engine nacelle or 
within a wing, the first term on the right side of equa- 
tion (13) vanishes. This fact makes the installation with- 
in a wing very attractive because, without the. frontal- 
area limitation, it is only necessary to seek a compromise 
between the second and the third terms. The best compro- 
mise results in an installation with a more favorable cool- 
ing efficiency than is possible when a separate nacelle in- 
stallation is used. 
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HADIATOES 
Apparatus for Heat Transfer 



The heat dissipation of several radiators was stud- 
ied. Circular tubes were used for the following radia- 
tors: 

(1) Length, 9 inches; tube diameter, 0.250 inch. 

(2) Length, 18 inches; tube diameter, 0*240 inch. 

(3) Length, 5 inches; tube diameter, 0.125 inch. 

In addition, a radiator was built' up of corrugated tubing 
(fig, 10) with a length of 1.97 feet and a mean diameter 
of 0.41 inch. for these four radiators, air was passed 
throtigh the tubes and water was directed perpendicular to 
them. The average inlet and outlet temperatures of both 
the air and the water were measured by means of thermo- 
couples. The mass flow of air was obtained either by a 
survey using a pitot tube and static tubes in the duct 
downstream from the radiator or, more commonly, by a ven- 
turi connecting the air at room temperature and atmos- 
pheric pressure with a large box, or expansion chamber, 
ahead of the radiator duct. The mass flow of water was 
obtained from a calibrated orifice in the water line. 
The pressure drop across the radiator was measured by 
pitot and static tubes connected to an alcohol, manometer. 
The apparatus for the heat-transfer investigation is shown 
in figure 11. 

In the results herein presented, the actual average 
value of the local heat-transfer coefficient from air to 
tube wall h a is used. In order to' obtain an accurate 
value of this coefficient, the water passages were made 
short and the water velocity was kept high so that the 
drop in temperature of the water AT W was very small. 

Several water speeds were used at each air speed, and the 
over-all heat-transfer coefficient h^ was plotted against 
AT w /AT a , as in figure 12, where each curve is for a con- 
stant mass flow of air. The curves were extrapolated to 
AT w /AT a =0. At this value, l/i^ = 0, where is the 

heat-transfer coefficient from liquid to tube wall. Since 
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JL = JL + JL- 
h t h a h l 

then h a = h t . The value of h t was computed in all 

cases using the logarithmic mean temperature difference 
from air to water. ; 

Heat-Transfer Data 

The formula for the heat-transfer coefficient given 
by McAdams (reference 3, p. 173) is for completely devel- 
oped turbulent flow. Figure 13 shows the heat-transfer 
data, plotted in nondimensional form, obtained from sever- 
al typical radiators compared with the McAdams formula. 
Except for the corrugated- tube radiator, the radiators 
tested showed a lower heat-transfer coefficient than indi- 
cated by the formula for turbulent flow. The same general 
formula can be used, however, for radiators at Reynolds 
Numbers above 10,000 by changing the coefficient from 
0.0269 to 0.0247. Similarly, Hartshorn (reference 2) ana- 
lyzed the data of Lorenz (reference 5) and found the- coef- 
ficient equal to 0.026. Care must be used in the follow- 
ing sections to apply the data and the equations to regions 
where they will be valid. 

The heat-transfer coefficient was found to be lower 
than that for turbulent flow given by McAdams and Harts- 
horn (fig. 13). This difference was to be expected inas- 
much as the friction loss was also less. (See fig. 4.) 
The explanation lies in the entrance conditions for the 
flow. At Reynolds Numbers in the usual operating range of 
a radiator, an appreciable length of tho tube is in the 
stabilization region where the flow is neither laminar nor 
turbulent. Thus the friction factor and the heat-transfer 
coefficient should fall between the curves for established 
laminar and turbulent flow. 

The empirical equations presented earlier based on a 
turbulent boundary layer will give both too high a heat 
transfer and too high a power consumption to push the air 
through the radiator. Almost the same optimum dimensions 
will be found for the radiator regardless of whether these 
equations or the present experimental results are used. 
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INSTALLATION WITHIN A WING OR AN ENGINE NACELLE 



Tlie high cost' of cooling commonly associated, with 
liquid cooling was due almost entirely to the radiator in- 
stallation. The radiator was located in the air stream 
without "benefit of cowling and with no means of control- 
ling the cooling. The drag was large, owing to the large 
frontal area, and the cooling power increased with the 
cuhe of the air speed. If the radiator gave adequate 
cooling in climb, too much cooling was available at cruis- 
ing. Radiators have been made retractable and have been 
mounted in the engine nacelle, in the wing, or in separate 
nacelles. Each of these arrangements has lowered the 
power required for cooling. The cost of cooling can be 
reduced to a fraction of that required in present instal- 
lations by a properly designed arrangement. 

When the radiator is installed within a wing or an 
engine nacelle, the power chargeable to the radiator is 
composed of two parts: the power actually used in cooling 
Pj, which is given by QAp/T}p; and the power used to 

carry and propel the weight of the radiator and water P^., 
which is given by e (C-p /C L ) W r 7 Q . The results of these 

considerations will be presented on the basis of optimum 
heat dissipation with respect to the power dissipated for 
various tube lengths and diameters. 

For the following case, several variables will be 
kept fixed in order to simplify the discussion. Later 
these variables will be changed one by one to ascertain 
their importance and effect. The design Ap is taken as 
25.6 pounds per square foot; the entrance and exit losses 
for each radiator are taken as 0.2 q^; the pump effi- 
ciency for the duct is taken as 100 percent; the ratio of 
open to total frontal area f is taken as 0.650; 
€ (C Dw /C L ) = 0.1; the fluid- is water, T w - T ia = 70° F. ; 

the tube-wall thickness is 0.005 inch; and the density, 
the specific heat, and the viscosity of the air are taken 
at 73.5° E. 

The necessary data are obtained from figures 4 and 13. 
Figure 4 is ttsed to give the variation of tube velocity 
with tube length for Ap = 25.6 lb. /sq.ft. and several 
tube diameters, as previously explained. Then, by means 
of equation (10) with Vj. from figure 6 and. the heat- 
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transfer data from figure 13 (equivalent to putting c 1 = 

0.0247 in equation (10)), curves of — f-- against 

ttD /4 

length for the assumed conditions are plotted in figure 14. 

Figure 14 differs from figure 8 in several respects. 
The data for figure 8 were obtained from the empirical 
equations for fully developed turbulent flow, whereas the 
data for figure 14 are the results of measurements under 
the actual flow conditions. In addition, the exit loss is 
included in the pressure losses used to obtain figure 14, 
but the exit loss was not considered for figure 8. 



The Selection of Optimum Radiator Dimensions 
The power expended for cooling is equal to QAp , or 



foot of open frontal area. 

The design Ap, 25.6 pounds per square foot, together 

with the velocities from figure 6, gives the power for 
cooling for any tube length and diameter. 

The power expended in carrying and propelling the 

weight of the complete radiator is a function of the flight 

speed, the aerodynamic characteristics of the airplane, 

and the type of radiator. That is, Pw = •€ (Cp /C^) ^ r *^ 0 » 

w 

€ (0-n /G T ) = 0.1; and V n = 200 and 300 miles per hour 
■"tut 

for the illustration. The remaining variable, ¥ r , may 
be approximately computed for round or hexagonal tubes as 
a function of L/D, of length, and of the free-area ratio. 
(See the appendix.) The method is found to check the 
weights of the radiators used by Harris and Oaygill (ref- 
erence 6) within 2 or 3 pounds per square foot of frontal 
area. This variation is about that encountered in manu- 
facture . 

If the heat dissipated is divided by the sum of the 
powers expended for cooling and weight propulsion for a 
particular tube diameter and length, figure 15 is obtained, 

showing the variation of •= — -— =— against length, for 



P D = QAp = 



(ttD 



/4) V t Ap 
ttD S /4 



= V f Ap foot-pounds per square 
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Figure 15 is the solution to the problem for the pre- 
viously mentioned conditions. Table II compares the op- 
timum values of these curves for an energy dissipation of 
500 horsepower. It is now of interest to investigate the 
effect of variations in these conditions. 



TABLE II 

A COMPARISON OP THE OPTIMUM RAD IATOR-DSS IGN 
SOLUTIONS OP FIGURE 15 



[Energy to he dissipated, 500 hp . j T w - Tj[ a = 70° P.; 
Ap, 25.6 lb. /sq.ft.; f, 0.650; c (C D /C L ), O.l] 



v • 

o 

(m.p.h. ) 


D 

(ft.) 


L 

(ft.) 


L 


H 

nl^74 
(hp./ sq. 
ft. open 
frontal 

area ) 


H 


Power used 

P D + P w 
(hp./ sq.f t . 

open 
frontal 
area ) 


Prontal 
area 


D 


P D +P ¥ 


required 
(sq.ft. ) 


200 


1/96 


0.58 


55.7 


111. 6 


11.92 


41.9 


6.89 




1/48 


1.105 


53.0 


112.6 


9.82 


50.9 


6.82 




1/32 


1.50 


48.0 


107.6 


8.31 


60.2 


7.14 


300 


. — 

1/9 6 


0.52 


50.0 


109.6 


9.83 


50.9 


7.01 




1/48 


.90 


43.2 


107.6 


7.96 


62.8 


7.14 




1/32 


1.30 


41.6 


102.5 


6.68 


74.8 


7.50 



The Choice of Ap 

The Ap chosen for the design is fundamental in de- 
termining the relative cooling efficiency of the radiator. 
In any particular problem it will be necessary to deter- 
mine beforehand the value of the pressure drop to be used 
in the radiator design. An upper limit is set by the 
climbing speed of the airplane. The lower limit is deter- 
mined by the amount of space available for the duct and 
the radiator. 
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The choice of the design Ap is made as follows; A 
particular tube diameter is taken as fixed. Then, for 
each Ap considered, the tube velocity, the heat trans- 
fer, the power used for cooling, and the power used for 
weight propulsion are calculated for a series of values' of 
L, The conditions assumed and the calculations are shown 
in table III. Table III is used to construct figure 16(a), 

H 

showing a plot of cooling efficiency. ■= — , ■ against 

P D + P W 

length for the various values of Ap. The optimum length 
and the cooling efficiency for each value of Ap from 
this figure are used to construct table IY, which compares 
the various values of Ap on the basis of the power and 
the frontal area required to dissipate 500 horsepower. 
Data for Ap = 6.4 lb./sq.. ft. are taken from figures 6 
and 14. Figure 15(b) is then drawn, the values from table 
17 being used. Prom figure 16(b) it is evident that the 
choice of Ap to result in an optimum cooling efficiency 
will be determined by the space limitations of the wing or 
the engine nacelle. Although a fixed pump efficiency has 
been used in these calculations, a decreased pump effi- 
ciency will probably be associated with large velocity re- 
ductions in the duct entrance section. A similar set of 
curves will result for any other tube diameter since the 
optimum values of figure 16(a) are primarily a function 
of L/D. Hence, the Ap chosen from considerations of 
one tube diameter will serve satisfactorily for the design 
pressure drop. 
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TABLE IV 

A COMPARISON OP TEE. OPTIMUM VALUES PROM PIGURE 16(a) 



[Energy to "be dissipated, 500 hp.; T w - Tj. a = 70° P.; 
f, 0.750; T|p, 100 percent; € (Cd w /C l ), 0.1; 
V 0 , 300'ffl.p.h. } D, 1/48 ft. 3 



Ap 

(lb. /sq.ft. ) 


L 

(ft. ) 


E 

P D + P W 


H. 
ttD^/4 
(hp.) 


Power 
required 
Pp + P ¥ 

(hp, ) 


Prontal area 
required 

(sq. ft.) 


6.4 


0.36 


10.70 - 


'38.5 


46.7 


13.00 


25.6 


.'98 


8.80 


111.2 


56.8 


5.99 


40.0 


1.25 


7.63 


144.0 


65.5 


4.63 


64.0 


1.74 


6.25 


187.5 


80.0 


3.55 


100.0 


2.00 


4.98 


237.2 


100.5 


2.81 



It is apparent from table IV that the cooling effi- 
ciency increases as the Ap decreases. These cooling 
efficiencies are, however, for radiators of -varying 
lengths. Accordingly, computations were made to ascer- 
tain the variation of cooling efficiency for a given radi- 
ator as the Ap changes over a wide range. Por this pur- 
pose, equations (6) and (8) were combined to give 

/ _ ^ 4h L j 

__!__ = gp Vt Cp (Tw _ TiJ \ ± _ e P ? t -Op € Dj 

where H is given in B.t.u./sec. 



V t , in f.p.s, 
p, in slugs/ cu. ft. 
g , in f t ./ sec. 3 



The data for V^ and h for each value of Ap are oh- 
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tained from figures 4 and 13, respectively. Figure 17 

H 

shows the cooling efficiency v~~T~%~ plotted against Ap 

for a standard-type radiator under the conditions given in 
the figure. It is seen that t'he cooling efficiency reaches 
a maximum at a A p of 13.5 pounds per square foot for the 
stated conditions, whereas the cooling efficiency for the 
optimum design at each Ap (fig. 16(h)) increased with 
decreasing Ap as long as the type of flow remained un- 
changed. 

The Variation in Pump Efficiency 

Figure 18 shows the effect of reducing the pump effi- 
ciency from 100 to 80 percent. In this case, the power 
expended for cooling, QAp/T|p f is increased; the efficien- 
cy of cooling is reduced; and the peak of the curve is 
moved in the direction of longer tube lengths. Harris and 
Secant (reference 7) present data on pump efficiency for 
wing ducts. 



Variation of the Ratio of Open to Frontal Area, f 

The heat dissipation per square foot of open frontal 
area is, of course, unaffected. The variation of f in- 
troduces two important changes, "both of which affect the 
weight of the radiator. An increase in f will decrease 
the water capacity of the radiator and the frontal area 
for any required quantity of heat dissipation. The maxi- 
mum value of f for a particular tuhe will he discussed 
on the oasis of a minimum desirahle water-passage width. 

Somewhat arhitrarily, the minimum water-passage width 
is set at 0.028 inch in order to assure open waterways and 
a moderate expenditure of power for pumping. The thick- 
ness of the tuhe metal is taken as 0.005 inch. The free- 
area ratio for these conditions can now he calculated from 
geometric considerations. (See tahle 7.) 
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TABLE 7 

THE EFFECT OF TUBE SHAPE AND SIZE UPOU THE FREE-AHEA RATIO 

[Wall thickness, 0.005 in.; minimum water-passage 

width, 0.028 in. ] 



Hadiat or 


Tube shape 


Inside diameter 
(ft.) 


Free-area 
ratio, f 


A 


Round 


1/96 


0.465 


3 


do . 


1/48 


.60 


C 


do . 


1/32 


.65 


D 


Hexagonal 


1/96 


. 5 9, 


E 


do . 


1/48 


.75 


F 


do . 


1/32 


• . 83 



The diameter for the hexagonal tube is the distance 
across the flats, equal to the hydraulic diameter. 

Figure 19 shows the variation of cooling efficiency 
with length for these radiators. Table 71 lists these 
radiators at the peak values. The corrugated- tube radia- 
tor is included for comparison. 

On the basis of this analysis, the hexagonal tube is 
shown to be superior to the round tube. Neither tube ap- 
pears to have any advantage from a consideration of ease 
of construction. 



The Effect of Variation in C D / C L 
The power used to carry the weight of- the radiator is 



duce exact 
in 7 0 . Figure 15 gives the plot of 



Thus a variation in Oj) /Oj, will pro- 

w 



duce exactly the same effect as a proportional variation 

H . 
— — — against 

P D + P W 

length for velocities of 200 and 300 miles per hour with 
€ (Oj /C^) constant at 0.1. The same figure can be used 

to show the variation of e (0^/0^) by keeping 7 0 = 200 
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miles per hour and giving € (Cd^/Cl) the value of 0.1 for 

the upper curve and of 0,15 for the lower curve or, alter- 
nately, by keeping V Q « 300 miles per hour and giving 
e (C-n /G L ) the value of 0.1 for the lower curve and of 

0.067 for the upper curve. Inasmuch as the value of 
e (°D / c l) is fixed "by the particular airplane considered, 
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it will affect the design considerations only as it affects 
the vaiue of the tube length for the maximum cooling effi- 
ciency. This fact is illustrated in figure 15, the curves 
being interpreted as suggested. • 



The Effect of Changes in the Air Constants Due to Heating 

The effect on the heat dissipate_d.- Consider equation 
(10) in the following form: 



E 



ttD 



Jl = P V t c p< T w - T ia> 



-4Ci -071? 



1 - e 



( P v t ) 



1 L 
T5T2 ~ 



(10) 



The change in c^ with ■ temperature is negligible. 

Consider the effect of a change in the value of the vis- 
cosity p.. A change in air temperature from 32° F. to 
122° F . will change p. from 0.0172 to 0.0193, about 2 per- 
cent in p. 0 * 2 . Thus, by taking an average value for p, 
the change would be about 1 percent in the exponent. The 
maximum value of H and the, L/D for the maximum value 
are not affected. For values of L/D below 35, however, 
the heat dissipation increases by about 1 percent. For 
values of L/D over 40, the effect is a fraction of 1 
percent and can be neglected. 

In equation (10), the density occurs only in the com- 
bination pY t . Therefore, if the mass flow of air through 
the radiator is kept constant, the variation in p due to 
heating will not affect the heat transfer because pT t 
will not change. This condition of' constant mass air flow 
will necessitate an increase in the pressure drop across 
the radiator with a corresponding change in the power re- 
quired for cooling. 



The eff ec t on the power for ..cooling . - For a heated 
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radiator, the pressure drop across the radiator is equal 
to the difference in total pressure (fig. 20): 

AP = P x + \ P^i 2 - (p a + | p s U 2 3) 

= P^Pg-f 2 p 2 V 2 2 - | P x v x S j (neglecting end losses) 



2 



■ p.-p. - ^ t' - 0 (") 



variation in p due to heating and pressure 
change along the tube is considered. The pressure effect 
"being present for a cold radiator, it will not "be included 
in the corrections. 

According to the momentum theorem, 

P x - P 3 = Ap f + (p a Y s s - Pi V^) 
where Apf is the drop in pressure due to skin friction 

Then / rp_ x , a 



A P = .A Pf + p x 7^ (| 2 - - l) (l 



+ an end loss 0.1 p Y 2 (16) 

2 a 

neglecting the pressure change along the tube. The change 
in Ap upon heating the radiator and keeping pV con- 
stant arises from the terms on the right-hand side of equa- 
tion (16). The increased velocity of the air in the tube 
will result in a larger pressure loss due to skin friction. 
The second term on the right represents the increase in 
pressure drop as a result of an addition of momentum to 
the air in passing through the radiator. A small portion 
of this added momentum is converted into total pressure, 
reducing the momentum pressure drop by a factor f 2 / 2 . 
Finally, the end loss slightly increases by reason of the 
higher value of the exit velocity. 

Lorenz (reference 5) finds that, for various values 
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of T w - T^ a , the value of Ap^/c^ remains constant and 
equal to that for the cold radiator. Then, since Ap f /q t 
is constant for the various radiator temperatures, 



heated 



mean 



mean 



using pY=pV=p V 

11 22 mean mean 



This relationship will 



give the increase in Ap^ for a given change in the den- 
sity of the air with sufficient accuracy for design pur- 
poses. Capon (reference 8) shows from theoretical consid- 
erations that, for values of L/D of about 40, the in- 
crease in the skin-friction term is about 0.3 the momentum 
term. Ho data are available for longer tube lengths. 

The magnitude of the effect may he computed for one 
of the radiators tested: Consider a radiator with L = 
5 inches, D = l/8 inch, L/D = 40, f = 0.6. Take the 
initial conditions at 73.5° 3T. with T x = 130 feet per sec- 
ond. Then 



130 f.p.s, 



p = 0.00231 
r i 

\ = 1.67X10 



~ 4 



and. taking AT a =27° F. 
using P 1 \=P a \ 



\ = 533° r. 

a bsolute 



Then 

qj_= 19.5 lb. /sq.ft. 

End losSj. = 3.9 lb. /sq.ft. 

Momentum change = PjVj. 3 

(it - x A l - y) 



T 2 =560° P. absolute 



p 2 =0.00220 



V =136.5 f.p.s. 
u a =1.84 X 10- 4 



q 2 = 20.5 Id. /sq.ft. 

End loss 2 =4.1 lb. /sq.ft. 

= 1.6 lb. /sq.ft. 



H x = 8,110 

-- = 1.10 (by fig. 4) 
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(Ap- ) •=■ 21.5 11). /sq.ft. ( A P^) =22.0 lb. /sq.ft. 

± heated 

Ap i total=25.4 lb./sq.f t. (cold) Ap g total=27.7 lb. /sq.ft. 

(heated) 

The calculated change in pressure is 9.0 percent. The ob- 
served increase is about 6 percent. 

The condition- of constant mass air flow, assumed to 
simplify the foregoing discussion, caused a change in the 
Ap required. From design considerations, the available 
pressure drop is the independent variable. It is neces- 
sary to start with this quantity to compute the mass air 
flow and, finally, the power required for cooling. 

From equation (16), it follows that the increase in 
pressure caused by heating the radiator is 

Figure 4 will provide the velocity of the air in the 
tube Y x for a given available Ap. Equation (17) will 

then give the increase in the pressure drop after correct- 
ing for the temperature rise expected. Figures similar 
to 6, 14, and 15 are then constructed, using the reduced 
Ap for a cold radiator. As previously explained, the mass 
flow of air calculated for the cold radiator at the re- 
duced Ap will remain as the mass air flow for the hot 
radiator at the given Ap. In the calculation of the Pp 
or V t Ap, the Ap used is the actual pressure drop and 

- — -. i n the computation of the cooling efficien- 
cy for figure 15, the power for cooling is added to the 
power for carrying the weight of the radiator. 



The Effect of a Change in Altitude upon the 
Hadiator Efficiency 

Meredith (reference 9) in an analysis of the effect 
of altitude on cooling efficiency found the density as a 
function of the temperature. This function varies with 
the coolant and the initial temperature difference. Any 
solution of the problem is based on the standard . atmo s- 
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phore (reference 10), which itself is standard. only under 
rarely realized specified atmospheric conditions. 

A general method of ascertaining the effect, of alti- 
tude upon the radiator performance for a supercharged en- 
gine is needed. lor the case of a constant duct-exit open- 
ing, it is possihlo to evaluate the decrease in heat dis- 
sipation under specific flight conditions for increasing 
altitudes. 



Equation (10) is used in the following form: 

o )' L 



-^j- = 2.75X10 4 p VV-T ia ) 



0.2 



— 0 . 0 0 9 6 o , _ , _ 

^pVt 3/ 2 

X ■* 6 i 



where 



XX 



is in B. t .u./hr ./sq. f t . of open frontal area 



it 33 / 4 

p, in slugs/cu.ft. 

7^., in f.p.s. 

|A, in lh./hr. f t . 

The pro Diem reduces to obtaining the variation of y. , 
p, T w - T i a > and - 7 £ with altitude; equation (10) is 

then u&ed for the heat dissipation. The value for the air 
density at altitude may he corrected for the adiahatic 
compression in the duct "before the radiator. The conver- 
sion of dynamic pressure to static pressure within the 
duct increases the air density and the temperature of the 
air. (The quantitative increase in the temperature of the 
air is later discussed in the section entitled The Jet- 
Propulsion Effect.) 

For the particular airplane, it is necessary to know 
the variation -of velocity with altitude and the velocity 
of the air in the tube for the design altitude and condi- 
tions, for a constant duct-exit opening, V^. is directly 

proportional to • T Q . Hence, the variation of and 
p?^. with altitude may he computed. 

The quantity T w - Ti a varies with altitude and the 
coolant used. By a comhination of the foregoing data, the 
heat dissipation may he ohtained for any altitude, as 
shown "ay the following example. 
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The radiator is designed for sea-level operation; 
D = 1/48 foot, L/D =50, Y t = 150 f.p.s. It is desired 

to determine the variation of heat dissipation with alti- 
tude at a constant exit setting for level flight at maxi- 
mum speed. The performance calculations give the varia- 
tion of flight velocity Y Q with altitude,- as shown in 
figure 21. Next, the variation of Y t with altitude is 
computed. In the computation of pY t at altitude, the 
value of p for the free .air stream is used. With water 
as coolant and T w taken at 27° F. "below the "boiling tem- 
perature at the particular altitude, T^ - T ia is easily 
obtained. Similarly, with glycol as coolant, T w is taken 
at 52° P. "below the "boiling point for a 97-percent glycol 
solution. The values for T ia are for the standard atmos- 
phere. The necessary calculations are shown in table YII. 
The change in heat dissipation with altitude is plotted in 
figure 22 for "both the water and the glycol radiators. 
The decrease in heat dissipation with altitude is not so 
serious above the rated altitude because, after that point, 
the required amount of heat dissipation decreases. 

"With a constant exit setting, the decrease in heat 
dissipation with increase in altitude is greater for a 
glycol radiator than for a water radiator. The reason for 
this difference is that the variation of air temperature 
with altitude makes a larger percentage change in the tem- 
perature difference between air and coolant for the water 
radiator than for the glycol radiator. 

The Ap required for a constant rate of heat dissi- 
pation was calculated by means of equation (12 ) and is 
shown in table 711. In addition, the required velocity 
of the air in the tube, the power for cooling, and the 
cooling efficiency for a constant rate of heat dissipa- 
tion are shown. An examination of these columns \?ill read- 
ily show the importance of designing the radiator for con- 
ditions approximating those to be encotintered in opera- 
tion. Por example, the power used for cooling will have 
doubled at 20,000 feet for the glycol radiator even with 
the assumption that the pump efficiency has remained con- 
stant. The increase in Ap- r /q,. 0 with altitude may cause 

difficulty in dissipating the required amount of heat af- 
ter a certain altitude. Table YII is not intended to give 
a comparison of the altitude performance characteristics 
for water and glycol radiators. The same design was chosen 
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for both coolants to facilitate computation of table VII, 
and neither design is optimum for the conditions assumed. 

The preceding calculations neglected the corrections 
to the inlet-air temperature and density for the adiabatic 
compression in the duct-entrance section. Where only rel- 
ative values for the heat dissipation at altitude are re- 
quired, the use of the uncorrected values introduces very 
little error. 

SHAPE OP THE INDIVIDUAL TUBE 



The tests on single tubes showed that square and rec- 
tangular tubes had the same frictional resistance in the 
turbulent region as round tubes of the same hydraulic diam- 
eter. It in safe to predict, then, that the heat dissipa- 
tion for the square, hexagonal, or rectangular tubes will 
also be the same as for the round tubes. Indeed, the heat- 
transfer equations given for single tubes may be appliod 
to the noncircular tubes by using their hydraulic diameters. 
The only difforencos that should appear with variously 
shaped tubes should bo in tho heat transfer from liquid to 
tube wall, the power required to force the liquid through 
the radiator, tho weight of tho radiator filled, the ond 
losses, and the free-area ratio. 

For example, radiators formed of hexagonal or rectan- 
gular tubes, which have uniform liquid passages, will be 
lighter on account of smaller liquid 'capacity, will require 
less frontal area, and may require more power to force the 
liquid through the radiator. Of these effects, the ones 
of greatest importance - liquid capacity and frontal area - 
have been considered under the section Installation within 
a Wing or an Engine Nacelle. 

Some radiators have been built with irregularly shaped 
air passages, such as the corrugated tubes shown in figure 
10. Such shapes are used in an attempt to produce greater 
mixing of the hot air near the tubo wall with the central 
core of air. Obviously, tho desired result can be easily 
accomplished in a variety of ways. 

The variation of the heat- transfer coefficient with 
Reynolds Number for the corrugated- tube radiator is shown 
in figure 13. In order to apply the formulas and to com- 
pare this radiator with the other types tested, each tube 
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is considered to have a square cross section with side a. 
Then a = 0.410 inch and L = 23.7 inches. In addition, 
figure 5 shows the isothermal friction factor for this 
radiator as a function of the Reynolds dumber. 

For a Ap of 25.6 pounds per square foot, the heat 
dissipation per unit open frontal area as a function of 
length is given by figure 14. It is evident that, for 
equal values of L/D or equal lengths, the corrugated- 
tube radiator is definitely inferior in heat-dissipating 
ability. 

The low value of the air velocity through the" tube 
reduces the Reynolds Number to such an extent that the ef- 
fect of the higher value of the heat-transfer coefficient 
for a given .Reynolds Number is nullified and- the heat 
dissipation is reduced. Table YI compares the cooling ef- 
ficiencies for a corrugated- tube radiator with several 
circular- and hexagonal-tube radiators under the same con- 
ditions.. 

CONDITIONS OF FLOW AROUND THE TUBES 



In the material in the section Single Tubes, the 
variables affecting friction loss and heat transfer were 
considered for the inside of the tube. The phenomena in- 
volved for the outside of the tube must also be consid- 
ered. The same fundamental principles rule in the trans- 
fer of heat from the liquid to the outside of the tube 
as from the tube to the fluid on the inside.' The liquid 
may flow perpendicular or parallel to the tubes. 

•The fact that the same phenomena of friction loss and 
heat transfer exist for the outside of the tubes as for 
the inside would appear to make conditions on the outside 
as much of a problem for study as conditions on the inside. 
The over-all heat- transfer coefficient h t from liquid 
to air is given by 




If either h a or h^ is small, the total thermal re- 
sistance l/h t is largely determined by that particular 
term. In most liquid-to-air radiators, i/h a is much 
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larger than l/h^ and l/h^ ' can "be neglected. It may 
easily happen that, with low liquid velocities, h^ will 
become snail enough so that the thermal resistance l/h^ 
will be important. 

In the foregoing radiator design, h a is taken equal 
to h t , and a method of computing the tube length and the 
number of tubes for each particular case is given. In ad- 
dition, a certain water-passage width is taken as optimum. 
Thus, for the chosen tube diameter, the frontal area of 
the radiator is fixed. Since no conditions have been pos- 
tulated concerning the arrangement of this required frontal 
area into any particular form, it follows that the width 
and the depth (fig. 23(a)) may be varied to fit a desired 
installation within the limitation that the product wd is 
constant (fig. 23(b)). Considerations of the water flow 
around the tubes, however, lead to the establishment of 
criterions for maximum permissible values for either w 
or d. 

The quantity of heat dissipated for a particular de- 
sign is determined by the over-all heat-transfer- coeffi- 
cient and the temperature difference between the two fluids. 
In order to maintain a constant heat dissipation for in- 
creasing values of w, it is necessary to maintain both 
h t and T w constant. Further consideration will show 

that the constancy of h^. necessitates the constancy of 
AT W . For example, suppose h^ constant as AT W decreases 
with a greater radiator width. The decrease in At w im- 
plies a decrease in the quantity of heat lost by the water, 
whereas the increase in the average water temperature T w 
increases the water-air temperature difference and hence 
the quantity of heat dissipated (= h t AT w _ a A c ) also in- 
creases. The further assumption is herein made that the 
volume of water flow is constant for the given design heat- 
transfer conditions. 

Ko investigations have been reported on the variation 
of the heat- transfer coefficient with Reynolds Number for 
liquids flowing perpendicular to tube banks of round or 
hexagonal tubes at the low Reynolds Numbers found in radi- 
ators-. The present heat-transfer data, when round tubes 
are used, indicate a Wusselt number of about 65 at a 
Reynolds Number of 880 based on the water velocity at the 
minimum passage width and the tube diameter for the length 
parameter. The Nusselt number may then be roughly assumed 
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to "be directly proportional to the R 0 yno Ids' Number . Thus, 
h-l may be calculated for any value of w. Since h a is 
known from the design solution, it is possible to estimate 
the effect upon h^ and, accordingly, upon the heat dis- 
sipated for any increase in w. 

For long liquid passages or greater radiator depths, 
h^ becomes increasingly greater than h a . As previously 
explained, the increase in h^ does not affect the over- 
all heat-transfer coefficient. The only limitation to be 
considered is the power expended in pumping the water 
through the radiator. 



The type of liquid enters the problem by its effect 
on the heat-transfer coefficient, that is, through the 
Nusselt, Prandtl, and Reynolds lumbers. Prom the rela- 
tions shown earlier, the specific heat, the thermal con- 
ductivity, the viscosity, and the density of the liquid 
are the most important variables. Brown and Barlow (ref- 
erence 11) show that water gives about 6 to 8 percent bet- 
ter heat transfer than glycol. Equal flow rates of water 
and glycol were used and no correction was made for the 
diminished air flow due to the higher temperatures at 
which the radiator was operated. Since the diminished air 
flow decreased h a by several percent and since higher 
liquid-flow rates are possible with glycol at equal pump 
speeds (reference 12) it appears that the dimensions of 
the radiator and the velocity of flow can be made more 
suitable to the glycol and can reduce the 6 percent as much 
as desired. 

The heat-transfer coefficient h^ from liquid to the 
tube wall will vary in the same manner as h a from air to 
tube wall. Thus 



where values of n = 0.3 and m = 0.5 are taken for flow 
below the critical Reynolds Number from McAdams (reference 
3, ch. VIII). Various experimenters have found values of 
n and m differing from 0.3 and 0.5 over a limited range. 
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(19) 
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The average values chosen will serve sufficiently well for 
the discussion. Then 



at 194° F. , 
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(20) 



Then 



. / N 0.3 / ^ x 0.5 / \ 0,7 

Plater (_1 ) ( _±_ ) ) =1.92 

glycol " V °' 62/ ^°- 45/ 

In the normal operation of a radiator, h-j, = 50 h a 
when water is the liquid. Then 

JL = JL + JL i__ + JL 

h t h l h a " 50 h a h a 

50 

H t = 5I h a 

Thus, if h^ "becomes several times as large, the effect on 
h t , which is the determining factor, is negligible. If 
h; is only half as large, h t = 25/26 h &t and the differ- 
ence is 2 percent in h t . TJaus , by the use of glycol, the 
over-all heat-transfer coefficient at the same volume of 
liquid flow is 2 percent less than for water. The availa-^ 
ble temperature difference T w - Tj a with glycol is sev- 
eral times that for water, depending on the inlet-air tem- 
perature, so that the net result is several times the to- 
tal heat transfer with the same air flow and liquid flow. 
The cooling efficiency may be still further increased, for 
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the total available duct space, being adequate for the wa- 
ter radiator at a certain Ap , will now permit the gly- 
col radiator to function at a lower Ap with the desired 
result. 

INSTALLATION IN A SEPARATE NACELLE 



In the consideration of the installation of the radi- 
ator in a separate nacelle, previous considerations of 
power used for cooling and weight support apply. It is 
only necessary to add the power lost owing to propelling 
the nacelle without any cooling air flow thro gh the na- 
celle. It follows from equation (13) and the previous 
discussion that, for a required amount of heat dissipation, 



H £ ttD 3 - ; Cp n q 0 7 0 Hr tt£ 
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550 
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550 f J H 
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(21) 



Subscript 6 refers to the free air stream and Pp 
and P^r refer to unit open frontal area. 

The Choice of Ap 

The cooling efficiency improved with a decrease in the 
pressure drop for a radiator mounted in a wing or an en- 
gine nacelle. A radiator mounted in a separate nacelle, * 
however, has an optimum Ap , which is fully as important 
as the optimum radiator dimensions. For any particular 
tube diameter and flight conditions, the comparison can he 
worked out in exactly the same manner that was previously 
followed in Installation within a Wing or an Engine Na- 
celle. The value of C-p^ for the particular nacelle con- 
sidered is used. For the present comparison, a well- 
designed cowling is assumed for which C D ^ is 0.12 (refer- 
ence 4). The power expended in overcoming the form drag 
per square foot of open frontal area for the conditions in 
table III with the air at standard density is, from equa- 
tion (21), 29,5 horsepower. This value is combined with 
the data of table III and figure 24 is constructed showing 
the variation of the cooling efficiency 
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H 



H 



550 f 



with length for several values of Ap, It appears that, 
for any given installation, one value of Ap will give a 
maximum cooling efficiency. 

The Choice of the Tube Dimensions 

Figure 25 shows the variation of the cooling effi- 
ciency H/P t with length for various tube sizes and shapes 
for three values for Cj)^, two of which allow for the in- 
terference effect "between the wing and the nacelle. 

Table Till shows the values of L/D for the peak 
cooling efficiencies from figure 25. Considerable latitude 
is permissible owing to the flat peaks of the curves. Ob- 
viously, since the interference effect is very favorable, 
the same L/D is nearly equally good for all the inter- 
ference factors computed. 

* » 

It may be noted that, when the interference effect is 
moderate or negligible, there appears to be an optimum 
diameter as well as an optimum length. 

For any required amount of heat dissipation, the power 
required for the radiator is obtained by dividing the re- 
quired heat dissipation by the cooling efficiency. 
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TABLE VIII 

HAD IAS OS LENGTHS REQUIRED EOR THE PEAK COOLING EFFICIENCIES 



E T w ~ T ia = 70 ° *-S A P» 25 - 5 Id. /sq.ft.; 

f, optimum for each tube shape and size; TL » 100 percent; 
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The actual cooling efficiencies presented in this re- 
port are somewhat misleading owing to the small value 
(70° F.) assumed for T w - Ti a . In actual practice, the 
value of T w - T ia will "be two or three times the design 
value, increasing the heat dissipation "by the same factor. 
Accordingly, the cooling efficiency will "be two to three 
times as large. A radiator designed from the values in 
this report giving a cooling efficiency of, say, 7 will 
actually give a cooling efficiency of 14 to 21; that is, 
for every horsepower of heat to he dissipated, l/l4 to 
1/21 horsepower will he required, or approximately 3.5 to 
2.4 percent of the brake horsepower will he used for cool- 
ing. The importance of a high-temperature coolant liquid 
is evident, for the cooling efficiency is directly propor- 
tional to the initial temperature difference. 
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The Cooling Efficiency for a Radiator in a Free Air Stream 

The radiator chosen for the comparison is made up of 
hexagonal tubes having a hydraulic diameter of 0.06 inch 
and a length of 3.74 inches. The heat-transfer and drag 
data are given by Parsons and Harper (reference 13) for - 
air speeds of 60 and 120 miles per hour. Extrapolation . 
was made to a flight speed of 200 miles per hour. The 
power used to transport the weight was calculated, as for 
the ducted radiator. It is found that, at 200 miles por 
hour, with T w - T ia = 70° F. and € (C D /0 L ) = 0.1, the 

cooling efficiency is about 1.7. Since the same radiator 
is designed to dissipate an equal or a greater amount of 
heat in climb at a much lower speed, it is assumed that 
the radiator is retractable so that this cooling efficiency 
will hold. For a nonretractable radiator, the cooling ef- 
ficiency would be much less in addition to the overcooling 
that would occur. This efficiency is about one-fourth that 
found in a separate nacelle and becomes rapidly less as 
tho velocity of the airplane increases. 

COOLING-POWER COMPARISON OF LIQUID- COOLED 
AND AIR-COOLED ENGINES 



It would appear that the liquid-cooled engine has sev- 
eral inherent advantages. The radiator can be made large, 
thus allowing cooling at a small Ap and, since the power 

is proportional to Ap , this power may be reduced to 
any desired value. It has been shown, however, that there 
is a limitation imposed by the power to transport the radi- 
ator. The liquid-cooled engine is often pictured as com- 
pletely enclosed in the wing or the nacelle with no form 
drag chargeable to the engine installation. 

The air-cooled engine, on the other hand, is described 
as located in a cowling with a large form drag. Now, if 
tho wing is large enough to house the liquid-cooled engine, 
it is also large enough to house the air-cooled engine. 
It has been shown (reference 14) that the interference ef- 
fect between wing and nacelle is favorable so that the ef- 
fective drag of the nacelle in front of a wing is only 40 
percent of the nacelle drag alone. When it is considered 
that, in these tests, the ratio of the nacelle diameter to 
wing thickness was large and that the conductivity of the 
engine was also large, it seems reasonable to assume that 
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tho interference on a modern engine with modern cowling 
located in a wing of thickness comparable with the engine 
diameter will he so favorable that the nacelle drag will 
almost disappear. It is probable that the 40 percent 
found in the tests of reference- 14 will be reduced to the 
order of 10 percent on a set-up of the nature assumed in 
this paper. 

It was found in reference 4 that the drag of a na- 
celle of 52-inch diameter at 100 miles per hour was 45 
pounds without cooling-air flow and that a streamline 
nose gave a drag of 32 pounds. It was further found that 
the engine could be cooled with an additional drag of 20 
pounds at this speed. The cooling drag will remain the 
same regardless of wing thickness. If the interference 
is favorable, as has been assumed, the form drag reduces 
to 4.5 pounds so that the total drag at 100 miles por hour 
is 24.5 pounds; that is, 5.3 horsepower for cooling and 
1.2 horsepower for form drag. 

With a well-designed cowling, the cooling power at 
higher speeds will not increase and may possibly decrease. 
At 300 miles per hour, there is still 5.3 horsepower for 
cooling and 27 times as much power, or 32 horsepower, for 
form drag. This form drag is independent of the engine 
power and depends only on the engine diameter. The cool- 
ing powor is proportional to the engine power, Tho values 
cited arc for. a 550-horsepower engine. Thus, a 2,000- 
horsepowor engine of tho same diameter at 100 miles per 
hour takes 19.3 horsepowor for cooling and 1.2 horsepower 
for form drag and, at 300 miles per hour, it takes 19.3 
horsepower for cooling and 32 horsepower for form drag. 
This total makes the power required for the engine instal- 
lation and cooling only' 2.5 percent of the engine power. 

The power to cool and to transport the cooling system 
is so small with either liquid- or air-cooled engines that 
it can hardly be a determining factor in the selection of 
one in preference to the other. 

If the engines are mounted in nacelles in front of 
comparatively thin wings, the liquid-cooled engine may 
have a streamline nose, thus having a . drag of 32 pounds 
instead of 45 pounds at 100 miles per hour. A separate 
nacelle. or installation for the radiator, however, will be 
required, which will add easily as much as this difference 
in drag. 
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'THE JET- PRO PULSION EFFECT 



The effect of heating an air stream after expansion 
in a duct and then contracting the duct hefore expulsion 
of the air is to convert some of the added heat energy 
into thrust. Both Meredith (reference 9) and Capon (ref- 
erence 8) have made estimates of the thrust so derived. 

Heither Meredith nor Capon considers the weight or 
the form drag of the radiator. If the weight 'and the form 
drag are considered, the net effect is a reduction in the 
power chargeable to the radiator. 

In "both of the foregoing cases and in the following 
discussion, losses occurring at the duct entrance and exit 
are not considered. The question of an efficient duct de- 
sign is extremely important in the consideration of the 
propulsive effect, for an inefficient duct will more than 
offset any possible thrust to he so derived. 

The path traversed hy the air through the duct shown 
in figure 26 is illustrated hy the path ABCD on the pv 
diagram in figure 27. The point-s A, B, C, and D corre- 
spond to the positions marked on figure 26. 

Path A3 represents the adiabatic compression of the 
air entering the duct. Path BC represents the heating 
at constant pressure along the duct. Path CD represents 
the adiabatic expansion as the duct contracts after the 
heating section. Finally, path DA, representing cooling 
at constant pressure, is equivalent to expelling the heated 
air at atmospheric pressure and taking in other air at the 
same pressure and atmospheric temperature. 

When no heat is added, the air traverses the cycle 
AB, BA, and no work is done. For the cycle ABCD, the 
work done is 
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(Po V C ~Pl T l) + Pi ^2 + Y ~~f (Pl V 2 -Ps^) 

+ P 0 ( T o - v 3 ) 



y 



(p. v + p v - p, v, - p V ) 



(22) 



Y 



- 1 



42 



Using pv = RT , 



W = — — — — (T + T -I ~ T ) 
Y - 1 2 o 1 3 

The heat added per unit mass is 



(23) 



(24) 



'The efficiency of the conversion of hoat energy into me- 
chanical energy is 
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which gives, using the adiabatic relation 
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and the relation R = c 
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(26) 



a result already given "by Meredith (reference 9). 



Host of the mechanical energy is recoverable as thrust 
energy. The ratio of the thrust energy to the total mechan- 
ical energy, per unit mass of air, is 



V (Y ~ V ) 
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For example, take Y Q = 300 miles per hour and V 3 = 328 

miles per hour; the useful thrust energy then represents 
95.5 percent of the mechanical energy obtained. 

In practice, the heat is added by a radiator that re- 
quires a certain pressure drop for the heat dissipation. 
The assumption of constant pressure heating must ho altered. 
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Path BC (fig. 27) is replaced by some path BC 1 , with a 
decrease in the area of the cycle indicating a reduced 
conversion of heat into mechanical energy. In order to 
ascertain the effect of radiator design upon the propul- 
sive effect, the heat will still ho considered added at 
constant pressure, this pressure being taken as the aver- 
age over the radiator. 

Then for V 0 = 300 miles per hour at sea level and 

standard atmospheric conditions, p 0 = 2,116 pounds per 

square foot and q Q = 230.4 pounds per square foot. The 

design Ap is taken as 25.6 or 64 pounds per square foot. 
The respective efficiencies are obtained, using equation 
(26) with V = 1.4: 



2.8 and 2.6 percent of the 
heat dissipated converted 
into mechanical energy- 
It has been previously shown that, for any installation, 
there is a Ap which will permit cooling at a .minimum 
power expenditure. In the radiator design, the only fac- 
tor that affects the energy recovery is this drop in pros- 
sure across the radiator; the lower the pressure drop, the 
greater the recovery. 3For the installation in a wing or 
an engine nacelle, the lowest possible value of the Ap 
is used for design. Hence a further decrease in the Ap 
to increase the energy recovery is impossible. For the 
installation in a separate nacelle, the optimum Ap is 
used for design. (See fig. 24.) Whereas a change in Ap 
from 64 to 25.6 pounds per square foot will change the en- 
ergy recovered from 26 to 28 horsepower for a 1,000-horso- 
powor dissipation, the cooling efficiency decreases by 15 
percent or about 12 to 15 horsepower for a good design. To 
attempt a compromise would reduce the energy recovery to 
a negligible amount with a real increase in the power ex- 
penditure. The conclusion follows that a good radiator 
design without consideration of the thrust will also be 
the best design when the thrust effect is taken into ac- 
count . 

As pointed out- by Meredith (reference 9), the adia- 
batic compression of the air entering the duct causes a 
rise in air temperature. The available temperature dif- 
ference T w - T ia will be less than the value using the 
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atmospheric temperature "by an amount depending on the air 
speed. From equation (25), assuming the air brought to 
rest, 



T ■ 



gives the upper limit for the temperature rise. Meredith 
gives the same result in the simple form 




CONCLUSIONS 



1. The design of a radiator is a function of the con- 
ditions under ^?hich the radiator is to operate. Hence it 
Is necessary to know "beforehand: . 



(a 
(b 

(c 

(d 
(e 

(f 
(g 
(h 

(i 



The quantity of heat to he dissipated. 

The available space for an internal installation, 
o r 

The drag coefficient of the nacelle when no cool- 
ing air flow is permitted, in the case of an 
external installation in a separate nacelle. 

The Ap available in climb. 

The type of liquid, water or glycol, and the value 
of T w - T^ a expected. 

The probable pump efficiency of the duct or the 
nacelle installation. 

The ratio of the lift and the drag coefficients 
of the wing. 

The total weight of the radiator and the addition- 
al structure for any radiator dimensions. 

The velocity of the airplane in level flight. 
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Such information is both necessary and, when combined with 
the ho at- transfer and the pressure-loss data of this re- 
port , is sufficient for an optimum radiator design. Par- 
ticular cases may add othor qualifications, such as maxi- 
mum possible radiator length, a maximum Ap for design to 
insure good cooling on the ground,' etc . 

2. The effect of each of the following factors on 
the cooling efficiency and the radiator dimensions was 
considered. 

(a) The pump efficiency:' A decrease in the pump 
efficiency of the duct will proportionally increase 
the power for cooling and decrease the cooling effi- 
,. ciency.. ... ' . .... . ,.. 

(b) The free-area • ratio : An. increase in the 
value of the free-area ratio will decrease the radi- 
ator weight and the total frontal area with a corre- 
sponding increase in the cooling efficiency. 

(c) The ratio e (C-n /Or): The power used to 

transport the radiator weight is directly proportion- 
al' to the value of * ($1)^/ ^L^ ' 

(d) The change in the air constants, c p , 

and p due to the heating up of the air: To main- 
tain a fixed rate of heat dissipation requires the 
condition of, a constant mass flow of , air , when the ra- 
diator is heated. This requirement in turn necessi- 
tates a larger Ap than is required by the same air 
flow for the isothermal flow., The 'design solution 
includes this effect by working out the! isothermal . 
case at a reduced Ap. The calculation of the ex- 
pected increase in Ap caused by heating the radia- 
tor is presented. 

(e) The change . in the heat dissipation with al- 
titude: For a supercharged engine in which the max- 
imum required heat dissipation is at the rated height, 
an increase in the air flow through the duct is re- 
quired if the radiator is designed for ground-level 
conditions.'* The radiator must be designed with a 
view to the rated altitude at which it is to operate 
because the proportion of the available AP required 
for cooling increases with altitude.. 
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(f) The width of the water passage: It is de- 
sirable to use the smallest possible water-passage 
width, conditioned upon the power used to pump the 
liquid through the radiator and the necessity of 
keeping the water passages open , in order to reduce 
tho radiator weight and to increase the free-area 
rat io . 

A variation in any of., the foregoing factors increas- 
ing or decreasing any part of the power chargeable to the 
radiator will cause a change in the optimum radiator di- 
mensions tending to bring the various opposing power con- 
siderations into balance again. For example, suppose the 
pump efficiency decreases. Then P-p increases and the 

optimum radiator design for the new conditions will have a 
longer length. The velocity' for any given Ap will be 
less and the value of Pp : will be decreased. 

3. The fundamental data presented on heat transfer 
and pressure losses for smooth.- tube- radiators' were ob- 
tained for the flow . conditions that actually exist in the 
tube section Containing the entrance.' A comparison with 
the established results for fully developed turbulent, 
flow showed that the heat transfer at a particular Ap is 
about the same for both cases with a larger power expen- 
diture required to force the air through the radiator for 
the actual flow. 

4 . An efficient radiator installation within a duct 
of nacelle where the quantity of flow' is controllable is 
superior to any installation where the- radiator" is exposed 
to the air stream, whether the radiator is retractable or 
shuttered. Relative cooling efficiencies for the instal- 
lations under similar assumed conditions were as follows: 



I&s t al lajb i on 



Eel at i ve ' 
cooling ;eff iciency 



Retractable radiator 



1.7 



Separate nacelle 



5.9 



llo wing-nacelle in- 
terference 



Separate nacelle 
Internal- duct 



S.O 



9.8 



40 percent C-r* 
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5. In the design of a radiator for an installation 
within a wing duct or an engine nacelle, it was found that 
a higher cooling efficiency is pbssible with a decrease 

in the : Ap across the radiator. HenCe , the first step in 
such a design is to select a value for the Ap that will 
require a radiator filling all the 'available frontal area 
in the duct space. This Ap will he the optimum and min- 
imum for the particular case. 

6. The optimum design Ap for an installation with- 
in a separate nacelle is the result of a compromise be- 
tween the powers used for cooling, weight propulsion, and 
form drag. This case approaches the design problem for a 
duct as the effective drag- coefficient of the nacelle de- 
creases and the power used to overcome the form drag be- 
comes of decreasing importance. On the other hand, as the 
form-drag power increases in relative magnitude, the fron- 
tal area must be decreased, requiring a larger Ap and a 
longer radiator, thus increasing the power for cooling and 
the weight propulsion. The best design is in tho region 
whore the opposing considerations arc of equal importance* 

7. The L/D- values for the solutions of a radiator- 
design problem, different tube diameters being used, fall 
very near each other. Hence, the length available for the 
radiator is a secondary factor since a reduction in the 
radiator design length merely calls for the use of a small- 
er tube diameter. In many cases, one tube diameter will 
give a slightly higher cooling efficiency than- is possible 
with any other diameter; this solution will be theoretic- 
ally the most desirable, although perhaps not practicable 
owing to the length necessary. 

8. A comparison was made between the circular tube 
and the hexagonal tube. It was found that: 

(a) The weight of the radiator required for a 
necessary heat dissipation is lower for a hexagonal 
tube because the water content of the radiator is 
less for equal minimum water-passage widths. 

(b) The free-area ratio is much higher for a 
hexagonal tube. Since both shapes have equal heat 
dissipating ability per square foot of open frontal 
area, the necessary frontal area will be less for a 
hexagonal tube. This fact may be used to advantage 
in the light of conclusions (5) and (6). 
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9. The higher temperature differences available "be- 
tween liquid and air, when glycol is substituted for wa- 
ter as the coolant, increase the total heat transfer sev- 
eral times,; depending on the inlet-air temperature. The 
decrease in the over-all heat-transfer coefficient when 
glycol is substituted for water is shown to "be several 
percent. The net result is an increase in heat dissipa- 
tion per unit frontal area of. the radiator. The conse- 
quent space saving for any required amount of heat dissi- 
pation may he used to decrease the power required to over- 
come the form drag and the jravrer "xcsed in propelling the 
weight or, alternately, the radiator size may be kept un- 
changed and the cooling -accomplished at a lower Ap with 
a saving in the power for cooling. In the design process, 
the balance between the various powers is automatically 
taken, into account by the selection of the optimum design. 

10. A comparison of the cooling efficiencies for the 
air-cooled and the liquid-<coo.led engines was made. It was 
shown that the power used for cooling and weight propul- 
sion Is of comparable magnitude for both cases and is too 
small in either case to be used as a determining factor in 
the selection of one engine over the, other. 

11. It was found that, at high. speeds, a mechanical- 
energy recovery is possible from the heat dissipated by 
the radiator in the duct.. In the ideal frictionle ss case, 
most of this mechanical energy is usefully converted into 
thrust energy. Good radiator design and energy recovery 
are compatible for an installation within a wing. Por the 
installation in a separate nacelle, the detrimental effect 
of changing the design to increase the thrust effect is 
relatively too great to permit any compromise. 



Langley Memorial Aeronautical Laboratory, 

latlonal Advisory Committee for Aeronautics, 
Langley -Field, 7a. , December 20, 1938. 
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APPENDIX • . 
■THE" CALCULATION OF RADIATOR ... WE IGETS 



The weight of a radiator mas' - be divided into the 
weight of the tubes , the water within the radiator, the 
headers (filled), and the solder used to "join the tubes at 
each end, . • 

1. The tube material is taken to he copper of spe- 
cific weight 555 pounds per .cubic foot and the tube-wall 
thickness is 0.005 inch. Then the weight of the tubes per 
unit open frontal area = 



555 X 2l001 x _TrDL_ 
12 - s 



ttD / 4 



(round tubes) 



555 X 0^05 x ZJjVL 

12» I — 

> vJ5_ jj s 



(hoxagonal tubes) 



• / « 4 !>. i'\ 



0.926 x ^ pound 



(either shape) 



2. The weight of the water in the radiator is equal 
to the density of the water multiplied by the volume re- 
maining after the volume of the tubes and metal is deduct- 
ed from the total volume.. If f is the free-area ratio, 
then l/f represents the total frontal area per' unit open 
frontal area. Thus, the weight of the water, per square 
foot of open area 



= 62.4 L 



/ 



1 
f 



ttD x 



- 1 - 



0^05 
. 12 



tt D /4 



(round tubes) 



= 62.4 

i 



(i 



1 - 



'O.s.005 
3D 



^ ' (round or hexagonal tubes) 



3, The data of references. 6 and 15 were analyzed to 
give the relationship between the weight of the headers 
•filled and the tube length. Figure. 28 show.s the data and 
the curve used to obtain casing weights for the present 
report. The radiators of Harris and Caygill (reference 6) 
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had a frontal area of 1 square foot. Then, per square 
foot of open frontal area.,, the casing weight from figure 
2 8 is multiplied by J~l/f on the : assumption that the 
casing size increases with the square root of the frontal 
area. A small error is introduced when the total weight 
per square foot of open frontal area is used to compute 
weights for radiators of increasing frontal areas, because 
the casing weight is now multiplied, by the ratio of the 
frontal areas instead of by the square root of this ratio. 

4. The weight of solder per square foot of open fron- 
tal area is taken as 1.5 pounds. 
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TABLE HI 

THE CALCULATIONS REQUIRED FOR A CHOICE OF Ap 
[Hexagonal tube; computations based on unit open frontal area; T w - T la = 70°F.; f, 0.750; 
•Hp, 100 percent; ((Ctj^/Cl), 0.1; V 0 , 300 m.p.h.; D, 1/48 ft.; T a , 73.5 0 F.J 



(ft. 



rrD a /4! P D 



ftp=25-6 lb./aq.ft.kp=40 lb./sq. ft. 



(lb.) (hp.) j(f .p.s4(hp. ): (hp.) tt .p.ej (hp.) 



"t *tD 8 A 



&p=64 lb. /sq.ft. 



(hp.) 



(f.p.s.) 



(hp.) 



APalOO lb. /sq.f t. 



(hp.) 



f .p. si (hp.) i hp. ) 



TrD»A 



0.75 

1.00 

1.25 
1.50 

2.00 



63. 4 5-07 
84. 0 6.72 



125 .4 
167.0 



g.314- 

10. 04 

13-35 



148-7 102.1 

129.0 '112.0 



113.2 
102.3 
87-9 



115-5 
117.9 
119.0 



6.92 
6.00 
5.26 
4.76 
it. 09 



187.0 123-9 



13.60 



240.0 



i6it.o 



|i3g.o 



145.0 ji44.o 
131.0 1M-7.0 
112.5 1^9-0 



11.93 210.5 
10.55 j 187.5 
9.52 170.0 
8.18 1W.0 



153.3 27-9 

i 

170.0 |2lt.5 

i 

180.0 Pi. 8 

i 

I85.O J19.S 

188. 3 !l7.0 

i 



304.0 
266.0 
239.0 

218.3 
188-5 



186-5 {55-3 
208.0 4-8. 4 

221.5143.5 

1 

231.0:39.7 

j # 

237-2 3^-3 



TABLE VI 

OPTIMUM COOLING EFFICIENCIES FOR VARIOUS RADIATORS 

[Energy to be dissipated, 250 hp.; tube-wall thickness, 0.005 in.; minimum water-passage width, 

0.028 in.; T w - T la = 70°F.; AP, 25.6 lb./sq. ft.; Ti p , 100 percent; e(C Dw /C L ), 0.1; V 0 , 200 m.p.h.J 



Radiator 


L 

(ft.) 


L/D 


Heat dissi- 
pated per 33. 
ft. of open 
frontal area 

(hp.) 


Requ lred 
frontal 
area 

(sq.ft.) 


S Vfelght of 
complete 
radiator 

(lb.) 


Power re- 
quired 
for cool- 
ing P D 

(hp.) 


Power re- 
quired for 
weight P w 

(hp.) 


Total 
power 
re ou lred 

(hp.) 


Cooling 
efficiency 
(neglecting 
frontal area) 
H 

*D + P W 


A 


O.52 


50.0 


110.0 


4-. 88 


215 


12.25 


11.50 


23-75 


10.50 


B 


1.04 


50.0 


111.5 




251 


13.15 


13.40 


26.55 


9.40 


C 


1.48 


47.4 


107.0 


3.60 


295 


14. 40 


15.75 


30.15 


8. 33 


D 


.56 


53.8 


111.2 


3. 81 


190 


11.70 


10.15 


21.85 


11.45 


E 


1.17 


56.2 


114.0 


2.92 


215 


12.05 


11.45 


23.50 


IO.70 


F 


I.65 


52.8 


IIO.5 


2-73 


240 


13.20 


12.80 


26.00 


9.61 


Corrugated 
f = 0.82 


I.30 


38.0 


72.7 


4.19 


291 


14.7 


15.5 


30.2 


8.28 



^Neglecting decrease in casing weight as the totnl radiator size increases. 











Water 






Altitude 
(ft.) 


q ° 

(lb. /sq.ft.) 


(lb. /nr. ft. ) 


P v t 

( .SlURS ) 


V T i s 
(°F.) 


H 

(with fixed 
exit setting! 

(B.t.u./hr. ) 


, Ap 

(with fixed 
exit setting) 


*Pr 
required 

for 
oonstant 
heat 
dissipation 




Ap 
*r 

"0 




V t 

re co ire d 


P D 




H 

P D 




ft. a sec. 


(id./ sq.ri;. ; 




lib. /sq.ft. ) 






(f.p.s.)- 




(hp.) 






a. 0 


120 


0.04-3 


0-357 


126 


6.10 x 10* 


42.0 




42.0 




O.35O 


150 




11.5 


20.9 


a. 


b. 5,000 


116 


.042 


.326 


135 


6.02 


40.6 




41.7 




• 359 


161 




12.2 


19.7 


b. 


c. 10,000 


111 


.041 


.294 


143 


5-83 


38.8 




42.5 




• 383 


175 




13-5 


17. g 


c. 


d. 15,000 


105 


.0395 


.264 


152 


5.59 


36.8 




43.g 




• 417 


192 




15.3 


15.7 


d. 


e. 20,000 


97 


.0385 


.23^ 


160 


5.31 


34.0 




44.9 




.463 


812 




17.3 


13.9 


e. 


f. 25,000 


S9 


.0375 


.206 


169 


5.02 


31.2 




46.1 




• 5lg 


234 




19.6 


12.3 


f. 


g. 30,000 


go 


.037 


•179 


178 


It (CI 


28.0 




47. g 




• 598 


262 




22. g 


10.5 


9- 


h 35,ooo 


71 


.036 


•153 


186 


4.25 


24. g 




51.1 




.720 


297 




27.6 


8-7 


h. 


i 1-0,000 


58 


.036 




178 


3.36 


20.3 




66.9 




1.15 




383 




46.6 


5-2 


1 • 














Olyooi 




















TABLE VII 






1 

! 


V T ia 
(°F.) 


H 

(with fixed 
exit setting) 

(B.t.u./hr. ) 


Ap 

(with fixed 
exit setting) 

(lb. /sq.ft. ) 


Ap r 

required 

for 
constant 

heat 
dissipation 

(lb. /sq.ft.) 


Ap 

r 

% 


Vt 

required 
(f.p.a.) 


P D 
(hp. 




H 
P D 




THE EFFECT OF ALTITUDE UPON RADIATOR PERFORMANCE 


a. i 
t 


231 


11.2 x 10 6 


1*6.5 




46.5 


0.388 




150 




12.7 




34.7 


a. 








b- ! 


243 


10. S5 


44.g 




47.7 




.411 




164- 




14.2 




31.0 


b. 


[Heat dissipation and Pp are In hp./sq. ft. of open 


•■! 

d. 1 


251 
258 


10.2 
9.50 


42. g 

40.5 




51.6 
56.3 




.465 
.536 




184 
20g 




17-3 
21.3 




25.4 
20.6 


c . 

d. 








e. i 


265 


8. 78 


37-4 




60. g 




.627 




242 




26.7 




16.5 


e. 


frontal area; T) p , 100 percent, assumed for all values 


f ! 

9- ; 


272 

279 


g.Og 
7-3* 


34.3 
30.9 




65.9 
71.9 




.740 
• 899 




267 
306 




32.0 
4o.o 




13. S 

11.0 


f. 
9- 








h. ; 


286 


6.54 


27.4 




go. 4 


1.13 




354 




51. g 




8.5* 


h. 


of the exit setting; D, 1/48 ft.; L/D, 50 J 


1 i 

1 


276 


5.20 


22.4 




104 


1.79 




454 




85.9 




5-1 


i. 
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FIGURE LEGEHBS 



(a) Static-tube connection. 

(b) Sample tu'be. 

Figure 1.- Tube for friction-loss study. 

Figure 2.- Velocity distributions for the l/2-inch-diameter 
circular tube. 

Figure 3.- Pressure drop along a smooth, tube with a stream- 
line entrance. 

Figure 4.- Summary of friction-loss data for single tubes. 



Eadia- 


bor s 


Single 


tube s 


Curve 


Length 
(in.) 


Tube 
diameter 
(in. ) 


L/33 


Curve 


l/d 


A 

B 
C 
G 


18 
9 
5 

23.7 


0.240 
.250 
.125 
.410 


75 

36 

40 
Corru- 
gated 
tube 


D 
E 
F 


75 
36 
40 



Figure 5.- Isothermal friction data for several radiators 
and comparison with single-tube data. 



Curve 


D 

(ft.) 


Ap 

(lb./sq. ft.) 


A 


1/32 


25.6 


B 


1/48 


25.6 


C 


1/96 


25.6 


D 


1/48 


6.4 



Figure 6,- Air velocity in the tube against length for sev- 
eral diameters. End losses included; T & , 73.5° F. 

Figure 7,- Analysis of single tubes.. 

Figure 8. Heat transfer against length for several tube 
diameters. T w - T ia = 70° F.; fully developed turbu- 
lent flow; Ap, 25.6 Ih./sq. ft. 
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figure 9.- Heat transfer against length for several values 
of Ap. T w - T ia = 70° F. ; D, 1/48 ft. 

Figure 10,- Throo-tubo section from the corrugatea-tubo 
radiator. 



A 
B 
C 
D 
E 
F 
G 
H 
I 
J 
K 
L 
M 
H 
0 
P 



copper-coil secondary shorted. 

water pump. 

"bypass. 

valve. 

expansion tank. " 
orifice for measuring water flow, 
static tuhe , water manometer, 
venturi for measuring air flow.' 
thermocouple installation, 
static tube , alcohol manometer, 
"baffle plate, 
duct wall, 
radiator, 
expansion "box. 

pitot tube , alcohol manometer, 
installations for nine thermocouples. 



Figure 11.- Apparatus for heat- transfer investigation, 
Figure 12.- Illustration of extrapolation to find h a from 



h, 



for one of the radiators. 



Curve 


Radiator 
. length 
(in. ) 


Tube shape 


Diameter 
(in. ) 


A 


23.7 


Corrugated 


0.410 


B 


5 


Round 


.125 


0 


18 


do . 


• .240 


D 


9 


do . 


.250 


E 


Recommended for fully devel- 




oped turoulent flow 






' McAdams 


preference 3, p. 173) 



Figure 13..- Correlation of the heat-transfer data, 
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Curve 


D 

(ft.) 


AP 

(1-b. /sq.ft. ) 


A 


1/32 


25.6 


B 


1/48 


25.6 


C 


1/96 


. 25.6 


D 


1/48 


6.4 


G 


Corru- 


25.6 




gated 






tube 





Figure 14. ~ Heat transfer against length for several tube 
diameters. T w - T ia ' = 70° F.; end losses included* air 

constants taken at 73.5° Fi • ■: 

Figure 15.- Cooling efficiency against length. T w - Tj_ a = 
70° F. ; Ap, 25.6 lb./sq. f t . ? f , 0.650; TL , 100 percent; 

e(C D /C L ) t 0.1; Y, , 200 m.p.h.; TL , 300 m.p.h. 

w 

(a) Cooling efficiency against length for several values 

of Ap. T w - T ia = 70° "F. ; f, 0.750; T| p , 100 per- 
cent; e(C Dw /C L ), 0.1; T 0 , 300 m.p.h.; D, l/48ft. 

(b) The choice of the design Ap. Data from .table IV. 
Figure 16.- The choice of . Ap for a duct installation. 

Figure 17.- The change in cooling efficiency with Ap . 

T w - T ia = 70 ° P "5 f > 6-750; % , 100 percent; e (C D /C L ), 
0.1; V , 300 m.p.h.; D, l/48 f % . ; L, 0.75 ft.; T a , 

73.5° F . • ' : . 

Figure 18.- Cooling efficiency against length for two pump 
efficiencies.. T w - T ia = 70° F. ; Ap , 25.6 lb. /sq.ft.; 
f, 0.650; c (C Dw /C L ), 0.1; V Q , 200 m.p.h.; D, l/48 ft. 

Figure 19.- Cooling efficiency against' length for various 
radiators with optimum free-area ratios. A, B, etc., 
as in table V; T w - T ia = 70° F.j Ap , 25.6 lb. /sq.ft.; 
Tip, 100 percent; e ( Cj^/ C L ) , 0.1;. T o ,200 m.p.h. 

Figure 20.- The fluid conditions for a radiator in a duct. 

Figure 21.- The effect of altitude on the flow through the 
radiator. Fixed duct-exit opening; supercharged engine. 
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Figure 22.- Illustration of the change in heat dissipa- 
tion with altitude for a fixed exit opening. 

(a) The cross-flow radiator. 

(b) Possible frontal shapes for a given value of wd. 
Figure 23.- The length of the path on the liquid side. 

Figure 24.- Effect of Ap on the cooling efficiency for 
an installation in a separate nacelle. T w - Tj_ a = 70°F. 
fi 0.750; T) p , 100 percent; € (C D /Cj,), 0.1; T Q , 300 

m.p.h.; P, l/48 ft.; 0 D , 0.12. 

( a ) Cj) , 0.12; no interference. 

(b) Oj) , 0.048; wing interference. 

(c) Ou^, °«°18; wing interference. 

Figure 25.- Cooling efficiency against length for an in- 
stallation in a separate nacelle. A, B, etc., as in 
table V; T w - T ia = 70° F. ; Ap , 25.6 lb. /sq.ft.; f, 

optimum for each tube shape and size; T) , 100 percent; 
€ (°D w / c l). O.lj V 0 , 200 m.p.h. 

Figure 26.- Frictionless duct, heating at constant pres- 
sure . " 

Figure 27.- The pressure-volume cycle for the flow through 
a duct. 

Figure 28.- The analysis of the data on radiator weights 
from references 6 and 15. 
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